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1. Introduction 

Jet-like events, first reported by Brueckner & Bartoe (1983), 
are characterised by non-Gaussian spectral line profiles. Dere 
et al. (1991) suggested that they are produced by bi-directional 
jets as a result of magnetic reconnection. To date, these jet- 
like events (often called explosive events) are mainly observed 
in spectral lines formed at transition region temperatures (Dere 
1994; Chae et al. 1998; Innes et al. 2001; Madjarska & Doyle 
2003), although observations of explosive events in chromo- 
spheric lines are also reported. For example, Madjarska & Doyle 
(2002) presented the temporal evolution of different temperature 
plasma using high cadence (10 s) observations obtained with the 
Solar Ultraviolet Measurement of Emitted Radiation (SUMER) 
spectrometer, and found a time delay in the response of the 

5 vi 933 A (2 x 10 5 K) line with respect to Ly 6 (2 x 10 4 K), with 
the Ly 6 line responding earlier. They concluded that the jet- 
like events may first appear at chromospheric temperatures. In 
follow-up work, Doyle et al. (2005) reported on a joint SUMER, 
Coronal Diagnostic Spectrometer (CDS) on board the Solar 
Heliospheric Observatory and TRACE imager study, confirm- 
ing the possibility that some jet-like events originate in the chro- 
mosphere. They further suggested that jet-like events could be 
divided into two types: one formed in the chromosphere and the 
other formed in the transition region. Some of the observed fea- 
tures are the result of spicules and/or macrospicules (Madjarska 

6 Doyle 2003; Madjarska et al. 2006), while others are the re- 
sult of high velocity flows in small loops (Teriaca et al. 2004). 
In more recent work, Madjarska et al. (2009) presented observa- 
tional data relating explosive events to a surge and demonstrated 
that the division of small-scale transient events into a number 
of different subgroups, for instance explosive events, blinkers, 
spicules, surges or just brightenings, is ambiguous, implying that 
the definition of a feature based only on either spectroscopic or 
imaging characteristics as well as insufficient spectral and spatial 
resolution can be incomplete. 

Several numerical models were developed to study jet-like 
events. Sarro et al. (1999) used a ID magnetic flux-tube model 
to simulate the temporal evolution of UV emission line pro- 
files, e.g. C iv 1548.2 A, in response to energy perturbations 
located below the transition-region. The maximum blue-shifts 
they obtained reach values of the order of 100 km s . Innes 
& Toth (1999) presented a 2D MHD study on jet-like events 
with different initial conditions, representative of different re- 
gions in which the reconnection occurs, e.g. the corona and chro- 
mosphere. Their conclusion was that high-velocity components 
in the profiles of lines formed around 10 5 K can be obtained in 
both cases, irrespective of the initial conditions. However, heat 
conduction was not included, and no brightening was found at 
the zero velocity position of the spectral line. In their model, the 
initial equilibrium state consists of two regions of oppositely di- 
rected magnetic field lines, with a narrow current sheet between 
the two regions. Magnetic reconnection at the current sheet is 
initiated by introducing localized anomalous resistivity. Roussev 
et al. (2001) carried out 2D MHD simulations, where jet-like 
events are formed during the process of magnetic reconnection. 
In their model, the initial magnetic field is parallel to the y-axis 
(vertical), and there is a thin current concentration formed along 
the y-axis. Magnetic reconnection is initiated by a localized in- 
crease of the magnetic diffusivity in the current concentration. 
Blue-shifts of the order of 100 km s were obtained. By using 
the same model, they further extended the work and performed 
simulations under different physical conditions (Roussev et al. 
2001a, 2001b). 



Yokoyama & Shibata (1995, 1996) performed 2D magnetic 
reconnection to study coronal X-ray jets using both oblique and 
vertical initial coronal magnetic fields. The temperature of the 
hot X-ray jets they obtained reach 3 times the coronal temper- 
ature. Moreno-Insertis et al. (2008) considered magnetic recon- 
nection triggered by flux emergence from below the photosphere 
using a 3D MHD model. Very strong X-ray jets with high tem- 
perature (3 x 10 7 K at the reconnection site) and high velocity 
(peak velocity 400 km s _1 ) were produced. In their model, the 
flux emergence is very strong (maximum field strength 3.8 kG), 
which plays a very important role in forming such strong jets. 
Isobe et al. (2008) focused on the process of magnetic flux emer- 
gence, and presented simulations of magnetic flux emergence 
driven by the upward convective motion. They found that small- 
scale horizontal magnetic fields could be produced even when 
the initial magnetic field is uniform and vertical. The horizontal 
magnetic fields emerging from the convection zone into the pho- 
tosphere undergoes magnetic reconnection with the background 
vertical field, which is a source of high-frequency MHD waves 
that may contribute to coronal heating or solar wind acceleration. 

Murray et al. (2009) presented another simulation of mag- 
netic flux emergence, where the long-term evolution of magnetic 
reconnection was initiated by flux emergence. A series of recon- 
nection reversals (or oscillatory reconnection) was reported. All 
the flux emergence studies mentioned above did not include heat 
conduction and radiative effects. The latter will reduce both the 
temperature and the velocity of the outflow jets. 

Litvinenko & Chae (2009) discussed magnetic reconnection 
at different heights in the solar atmosphere, and found that the 
temperature and speed of the outflow jets vary by several or- 
ders. Their study was based on an extended Sweet-Parker model 
(Parker 1957, 1963), assuming that the inflowing magnetic en- 
ergy is completely converted in the current sheet into the thermal 
and kinetic energies of the outflowing plasma. 

In the present study, we use a 2. 5 -dimensional resistive MHD 
model in Cartesian coordinates to investigate magnetic recon- 
nection in the low atmosphere, e.g. chromosphere. We discuss 
its implications for jet-like events. Here, we use various phys- 
ical environments representing the lower atmosphere with dif- 
ferent magnetic field strengths and densities. In our model, the 
magnetic reconnection is triggered by newly emerging magnetic 
fluxes from below the photosphere which reconnect with the pre- 
existing background magnetic field lines in the chromosphere. 

The physical model is described in Section 2, while the nu- 
merical results are shown in Section 3. No background heating is 
included in this model, the implications of which are discussed 
in Sect. 4. 



2. Physical Model and Numerical Methods 

2.1. Basic Equations 

For 2.5-dimensional MHD studies in Cartesian coordinates 
(x, y, z), one may introduce a magnetic flux function t/r(f , x, y) to 
express the magnetic field by 



B = v x (if/z) + B z z. 



(1) 



The 2.5-D resistive MHD equations are in the following non- 
dimensional form 

% + v ■ (pv) = 0, (2) 
ot 

c h. + v . vv + - v p - -j x B - g = 0, (3) 

ot p p 
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where p is the mass density, v is the flow velocity, T is the tem- 
perature, p = pT is the gas pressure, g is the gravitational accel- 
eration, y(= 5/3) is the adiabatic index, r\ is the dimensionless 
magnetic diffusivity which is inverse to the magnetic Reynolds 
number, Bq is the characteristic ratio of the gas pressure to the 
magnetic pressure, j = V x B is the electric current density, 
Q = v ■ [r 5/2 (B ■ VT)B/B 2 ] is the field-aligned heat conduc- 
tion function, and L r is the radiative losses. Here we assume that 
the plasma is fully ionized, therefore, the dimensionless form of 
L, is explicitly given by 



L r = p z A(T), 



(7) 



where A(T) is the radiative loss function, of which the expres- 
sion used is that given by Rosner et al. (1978). Various physical 
environments representing the solar atmosphere with different 
characteristic values are examined. The different cases can be 
divided into two groups according to the characteristic value of 
Bq. The two groups are denoted by Exp. A and Exp. B, corre- 
sponding to /?o = 0.33 and 0.033, respectively. For both groups, 
the following characteristic values are kept invariant and taken 
as basic units: po = 3.33 x 10~ 10 kg rrT 3 for the mass density, 
Tq = 10 4 K for the temperature, Lq = 500 km for the length, and 
vo = ~\/RTq = 12.8 km s _1 for the velocity, where R is the gas 
constant. In terms of the expression Bq = 2fipoRTu/ B 2 where p 
is the vacuum magnetic permeability, we can derive the charac- 
teristic values of the magnetic field strength Bo = 6 G for Exp. A 
and 20 G for Exp. B. The dimensionless coefficient of the heat 
conduction and radiative losses are given by 



Ci 



(7 - l)K Tl /2 /(poL vl) and C 2 



(y 



l)L po/(vlml), 



respectively, where kq — 10" 11 W irT 1 K~ 7 / 2 is the classical con- 
ductivity coefficient, and m p is the proton mass. In the present 
study the x— and z-axes are in the horizontal plane, and y-axis 
represents the height of the solar atmosphere. Therefore, the di- 
mensionless form of the gravitational acceleration is expressed 
as 

g = ~g% 
where g = 0.271 km s _2 Lo/Vq. 



the chromosphere to the upper transition region, the profile of 
the initial temperature T is given by 

r(0, x,y) = = 1 + Ar{l + tanh [e(y - y,)]} 

where 2 AT - 100 defines the relative temperature jump across 
the transition region, e = 0.6 describes the steepness of the tem- 
perature profile and y s — 5 is the position of the middle tran- 
sition region. The temperature is equal to 1 at the bottom and 
is 101 at the top. Then, the mass density is calculated from the 
hydrostatic equilibrium equation, Vp = -pgy. In order to obtain 
a unique solution for mass density, we specify the density value 
at the bottom as a boundary condition. It is easy to deduce that 
both the temperature and the mass density of the initial state are 
dependent only on the y-coordinate, and are uniform along the 
x-coordinate. Therefore, the initial mass density at the bottom 
can be expressed by one parameter, p/,o. Fig. 1 shows one exam- 
ple for the initial distribution of the temperature and the density 
along the y-coordinate when p/,o = 3. 
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2.2. Initial State 

The initial magnetic field is a linear force-free field taken to be 
in the following form 



((4o)/tt) cos(7tjc/2w) - 4oj/tt) ^j2ppoRTo/Bo, 
(-2\x\ - 4co/n + 2u) ^2 m) RT //3 , 



B- 



2 cos(7rx/2w) ^2ppoRTo/Bo, 
0, 



ifW 
if W 



ifW 
ifW 

< 0} 

> OJ, 



< OJ 
> OJ, 

(8) 
(9) 



where oj is the half width of the force-free region (current sheet) 
and is set to 0. 1 . The unit of magnetic field strength is taken as 
Bq — sj2ppoRTo/Po. The lower /3q is, the stronger the magnetic 
field strength. In order to represent the solar atmosphere from 



Fig. 1. The distribution of the initial temperature (T, solid line) 
and mass density (p, dashed line) with height (y). The density 
distribution corresponds to the case where pt,o = 3. 



In this study, two initial states with different values of pbo are 
considered for each of Exp. A and B. In total, four cases which 
represent four physical environments are discussed here, all of 
which are listed below: 

- Exp. Al: Bq = 0.33, pw) = 3 (weak magnetic field and high 
mass density) 

- Exp. A2: Bq = 0.33, pbQ = 0.1 (weak magnetic field and low 
mass density) 
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- Exp. Bl: y6o = 0.033, pio = 3 (strong magnetic field and high 
mass density) 

- Exp. B2: y8o = 0.033, pho = 0.1 (strong magnetic field and 
low mass density) 

In all the cases, the dimensionless form of the initial magnetic 
field is the same, and the change of the magnetic field strength is 
only due to the change of Pq. Moreover, the dimensionless form 
of the magnetic flux emergence is the same in all the cases, as 
well as the dimensionless form of the magnetic diffusivity, which 
will be discussed in the following section. 

We should note that most numerical studies consider an ar- 
tificial hearting source, however we do not consider it here. The 
main reason is that we do not know where that artificial heating 
term comes from, thus, in our model, the initial solar background 
atmosphere is not in thermal equilibrium. In order to reduce the 
influence of the flows that are driven by the re-distribution of 
the heat and radiative cooling from the initial background at- 
mosphere, the calculations in our paper are carried out over a 
shorter interval (about 80 s), compared to the time scale for heat 
conduction and cooling from the initial configuration. 

2.3. New Magnetic Flux Emergence 

Here, we do not discuss the specific physical mechanism for the 
flux emergence which is beyond the scope of this paper, although 
recent observations with the Hinode Solar Optical Telescope 
(SOT; Tsuneta et al. 2008) reveal that the emergence of granular- 
scale magnetic flux is ubiquitous on the solar surface (Centeno 
et al. 2007; Ishikawa et al. 2008; Isobe et al. 2008; Okamoto 
et al. 2008). Instead, what we are interested in is the physi- 
cal response of the chromosphere after the magnetic flux has 
emerged from below the photosphere and reached the upper at- 
mosphere, the chromosphere and transition region. The emer- 
gence of new magnetic flux is implemented numerically by 
changing the boundary conditions of the magnetic field at the 
bottom, which is described below in detail. First, we specify the 
flux emergence region that is taken to be 1 < x < 4. Then, we 
change the magnetic flux function at y = in the emergence re- 
gion to a new value, ij/ e , according to a function of x and time, 
expressed by 

f e = (/%(*) + a\{*jr m {x) - lAin(l)) X (t]/ in (x) - <Ain(4))|f/f e , 

1 < x < 4, < t < t e , (10) 

where tfr m is the initial magnetic flux function at the base, t e is 
the duration for the flux emergence, and a controls the magni- 
tude and orientation of the emerging flux. When t is greater than 
t e , the newly specified flux distribution at the base remains in- 
variant in the emergence region. The flux distribution outside 
the emergence region is fixed to be the same as that of the ini- 
tial field, see equation (8). In all numerical examples discussed 
below, t e is taken to be 80 s and a is -1 .2. 

As new flux emerges, a current sheet will be formed at the 
interface between the newly emerging and pre-existing magnetic 
fluxes. Then, anomalous resistivity is introduced to initiate mag- 
netic reconnection. The distribution of the resistivity is localised, 
and taken to be of the following form 

( 77 m min(l,;7;' c ), if |*| > l,y < 3, 
' 1 0, elsewhere, ' 

where r] m = 0.002, and j c = 10, is the critical value of the current 
density (J). No magnetic reconnection is considered at the cur- 
rent sheet of the initial background magnetic field within \x\ < ai, 
see equation (11). 



There are thus two current sheets in our model. One is at 
the interface between the newly emerging magnetic flux and the 
pre-existing background magnetic flux and the other is along y at 
around x = which belongs to the background magnetic field. 
This is because we only want to study the magnetic reconnec- 
tion occurring on the former one. Here, we introduce a localized 
resistivity (at |jc| > l,y < 3). The current is almost outside this 
region (|x| > l,y < 3), except for that around x — 0. 

The initial state used here contains a current sheet around 
x — 0, which is not the same as the one used by Yohokyama & 
Shibata (1995, 1996) (where the initial magnetic field is a uni- 
form vertical one). The reason for using non-uniform magnetic 
field is that it is easier for us to compare the results in this pa- 
per with those where reconnections at both current sheets are 
considered. The two current sheets represent different magnetic 
topology: one is formed between open vertical field lines while 
the other is between open and closed ones. The reconnection be- 
tween the newly emerging and pre-existing magnetic fluxes will 
have influence on the reconnection between the open vertical 
fluxes and vice versa. In this case, the polarity of the emerging 
flux will be an important factor. However, in this paper, based 
on the fact that no magnetic diffusion is considered at the initial 
current sheet around x = 0, the influence of the initial current 
sheet on the results can be regarded as the influence of different 
left-boundary conditions if compared with the uniform vertical 
field case. As the reconnection occurs on the right-hand side of 
the emerging flux arcade, far away from the initial current sheet 
around x = 0, there is no significant differences in the results if 
a uniform initial field is used instead. 

2.4. Computational Domain 

The dimensionless size of the computational domain is -5 < 
x < 5 and < y < 10. Because of symmetry, we only implement 
calculation in the right half region (0<x<5,0<y< 10). As 
for the boundary conditions, we fix all quantities at the base, and 
treat the top (y = 10) and the right-hand side (x = 5) as open 
boundaries. Furthermore, we use symmetrical conditions for the 
left boundary (x = 0). A multistep implicit scheme (Hu 1989) is 
used to solve equations (2)-(6). 

We adopt 399x400 grid points, which is based on the consid- 
eration of the balance between resolution and computation time. 
We can also adopt more grid points so as to get better resolution, 
but the computation time will become very long then, especially 
for the case in which strong magnetic field strength and low mass 
density are presented. Uniform meshes are adopted in both the 
x— and y— directions. 

3. Numerical Results 

We present here magnetic reconnection which occurs in the low 
atmosphere and is triggered by newly emerging magnetic flux. 
Four cases representing different physical environments are dis- 
cussed. According to the strength of the magnetic field, the re- 
sults are divided into two groups: Exp. A and B. For each group, 
two initial states with different mass density at the bottom, p^o, 
are considered. 

3.1. Exp. A (fio = 0.33) 

First, let us discuss the case where p^o = 3 (i.e. electron density 
N e = 6 x 10 11 cm -3 ), which means the mass density at the bot- 
tom is large, corresponding to Exp. Al. The results are shown 
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in Fig. 2, where the evolution of the magnetic field, tempera- 
ture and velocity are displayed. At the beginning of the exper- 
iment, downward flows are seen, especially in the middle tem- 
perature region (see Fig. 2b), which is caused by heat conduc- 
tion and radiative cooling. As the new flux emerges, an electric 
current sheet is formed at the right interface between the newly 
emerging and pre-existing magnetic fluxes. Then magnetic re- 
connection occurs at the current sheet producing bi-directional 
outflows: upward and downward. The outflow jet in each direc- 
tion will be accelerated by magnetic tension forces and shocks 
along the boundary of the jet (Petschek 1964), and reaches its 
maximum velocity (vj et ) somewhere in front of the jet. For ex- 
ample, at t — 32 s, the upward flow jet reaches its maximum ve- 
locity, Vjet = 40 km s _1 , at x — 3.2, y = 1.1, whereas the diffusion 
region is located at a lower altitude, around x = 4.0, y = 0.3. As 
the magnetic reconnection continues, the upward flow jet propa- 
gating towards higher altitudes will interact with the downward 
flow of the background field, forming a shock at the interface 
between them. The shock will impede the downward flow, so 
that the downward flow first becomes weaker and then even re- 
verses its direction, see Fig. 2. During the process of magnetic 
reconnection, both the magnitude and orientation of Vj et change 
dynamically. But, no matter how Vj e t changes, it is far smaller 
than that of the jet-like events observed by SUMER. In spite of 
this, it is comparable with the characteristic local Alfven speed 
around the diffusion region, as predicted by the fast magnetic re- 
connection model (Petschek 1964; Priest & Forbes 1986). In the 
following sections, the velocity Vj et we give is the maximum one 
during the time from - 80 s. Also, the velocity in the line-of- 
sight relative to Vj e t will be presented. In this case, the downward 
flow jet (as a result of magnetic reconnection) can not be clearly 
shown, because the X-point is very low in our model. Its height 
is 0.2 above the bottom at t = 32 s, and is 0.8 at t = 80 s. 
The height is so small that only a small part of the downward 
flow jet is displayed in our computation domain. Moreover, as 
mentioned above, the outflow jet will reach its maximum ve- 
locity somewhere in front of the diffusion region, which means 
the maximum velocity of the downflow jet may be located below 
the bottom of our computational domain. The downward flow jet 



that can be displayed in our model reaches 12 km s at t — 80 s. 
In the following discussions, we will focus on the upflow jet. 

Fig. 2 is drawn in such a large scale that the temperature evo- 
lution around the magnetic reconnection region has a very poor 
resolution. In order to better illustrate the temperature evolution 
in this region, we plot only a small part (1 < x < 5, < y < 1 .6) 
of the computation domain in Fig. 3 (and in successive plots), 
where the relative changes in the temperature, (T - T m )/T m , with 
respect to the initial temperature, T m , are shown. The relative 
changes in mass density, (p - pi n )/p\ n , with respect to the initial 
density, p m , are shown in Fig. 4 at the same times and in the same 
region as in Fig. 3. 



t = 16.0 s t = 32.0 s 




Fig. 3. Relative changes in the temperature, (T - T\ n )ITi n for 
Exp. Al, with respect to the initial temperature (Tin), are shown 
at four times in the region 1 < x < 5, < y < 1.6. 
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t = 16.0 s 



t = 32.0 s 





t = 64.0 s 



t = 80.0 s 




Fig. 4. Relative changes in mass density, {p-p m )lp m for Exp. Al, 
with respect to the initial density (p m ), are shown at the same 
times and in the same region as shown in Fig. 3. 



Peak structures are clearly seen on each rim of the flux emer- 
gence region for both the relative changes in densities and tem- 
peratures. As the magnetic reconnection occurs, the plasma in 
the diffusion region is heated, and the hot plasma is pulled out by 
the magnetic tension force in either direction along the current 
concentration, resulting in the peak structure on the right-hand 
side of the flux emergence region in both Fig. 3 and 4, compared 
with those in the inflow region. Along the peak structure, the 
temperature reaches its maximum in the diffusion region, while 
the density is minimum in the diffusion region. After the mag- 
netic reconnection, part of the reconnected magnetic field lines 
extend outward towards x = (see Fig. 2), resulting in a negative 
horizontal velocity component of the upflowjet. Accordingly, on 
the left-hand side of the flux emergence region, plasma accumu- 
lates and a strong shock is created which causes the formation 
of the another peak structure in the relative changes in density 
plot; the peak structure in the relative changes in temperature 
plot is due to the heating effect of the shock. In the following, 
the temperatures at two locations are analysed: one in the diffu- 
sion region and the other at the location of Vj e t- They are denoted 
by T m . dx and 7j et , respectively. 

Even though the plasma temperature along the current con- 
centration is heated through magnetic reconnection, the relative 
changes in temperature are negative all over the region shown in 
Fig. 3. Both the temperatures r max and 7j et are very low, around 
10 4 K. In the following, we make quantitative analysis of how 
this happens. 

From equations (6) and (7), we get 



dT 2(y- I);? , 

- f - C 2 pA(T), 

dt P B J W 



(12) 



where the contribution of the Joule dissipation and radiative 
losses are considered. For simplicity, the first term on the right- 
hand side of Eq. (12) is denoted by E\, and the second term by 
Ei. Because both the temperature and its gradient are very small 
at low temperature regions, the heat conduction term (~ 10 -5 ) 
is neglected here. Note that the Joule dissipation and radiative 
losses play opposite roles in the evolution of temperature, one 



heats the plasma and the other cools it. From equation (12), we 
find that E\ is inversely proportional to the density, whereas E% 
is proportional to the density. In the diffusion region where the 
density is small, E\ is ~ 80, and Eo is ~ 50. The net effect is that 
the plasma in the diffusion region is heated. However, the heated 
plasma is pulled out of the diffusion region and the cold plasma 
outside the diffusion region is pushed in, so that the plasma in the 
diffusion region cools. Outside the diffusion region, the plasma 
density along the current sheet increases by several times with 
respect to the initial density (see Fig. 4), so that E% increases 
to ~ 200 at maximum, resulting in a fast temperature decrease. 
If the plasma temperature is below 10 4 K, radiative losses are 
reduced to zero, so that the relative changes in the plasma tem- 
perature is flat along the outflow region (located along the "high" 
temperature ridge on the right-hand side of the flux emergence 
region seen in the bottom right panel of Fig. 3). Furthermore, 
outside the diffusion region, the outflow plasma will expand, 
which will cool the plasma. The contribution of the plasma ex- 
pansion to plasma cooling is (y - l)T V -v ~ 40. In the diffusion 
region, the plasma expansion will also cool the plasma there, 
however its contribution is very small, only about 4. A large 
plasma cooling region is found around x — 2 - 4 almost reach- 
ing y = 2 (see Fig. 3), located in the inflow region. This is be- 
cause of the plasma expansion in the inflow region, of which the 
contribution to plasma cooling reaches 40. Even though the con- 
tributions of plasma expansion to plasma cooling are of almost 
the same order in the inflow region and the outflow region, the 
plasma temperature decreases more strongly in the former than 
in the latter. This is caused by the fact that there is no hot plasma 
flowing into the inflow region, whereas hot plasma heated in the 
diffusion region is expelled out into the outflow region. 

It seems that small-scale jets can not originate in a physical 
environment like Exp. Al where high mass density and weak 
magnetic strength are presented, as both the velocity and the 
temperature of the outflow jets are very low. What will be the re- 
sult if the physical environment with the same magnetic strength 
but a lower mass density is considered? The analysis shows that, 
if the density is lower, more plasma can be heated (because both 
the radiative losses and the ratio of the radiative losses to Joule 
dissipation is proportional to p 2 ). 

Fig. 5 shows the evolution of the temperature for the case 
of Exp. A2 where Bo = 0.33 and p/,o =0.1 (this corresponds 
to magnetic field strength B — 6 G, electron density N e = 
2 x 10 10 cm -3 ). Significant relative changes in temperature are 
clearly seen in this case and the plasma in the diffusion region is 
heated to the temperature of T m - dx = 2.8 x 10 5 K. The maximum 
velocity of the outflow jet caused by magnetic reconnection and 
its temperature are also analysed for Exp. A2. We find that Vj e t 
reaches up to 85 km s and 7j et is about 2 x 10 5 K. Both are 
larger than that obtained in Exp. Al. The direction of vj et is 32° 
with respect to the line of sight. 



3.2. Exp. B(fio = 0.033) 

In this section, the magnetic reconnection takes place in a phys- 
ical environment with a strong magnetic field. Fig. 6 shows the 
evolution of the temperatures for the case of Exp. Bl where 
Bq = 0.033 and pbo = 3 which corresponds to B = 20 G and 
N e — 6x 10 11 cm -3 . This represents a physical environment with 
a strong magnetic field and a high mass density. Fig. 7 shows the 
relative changes in density. Peak structures are present in the 
relative changes in both the temperature and density on each rim 
of the flux emergence region. In the diffusion region, T max is 
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Fig. 5. Relative changes in the temperature, (T - T m )IT m for 
Exp. A2, with respect to the initial temperature {T m ), are shown 
at four times in the region 1 <x<5,0<y< 1.6. 



Fig. 7. Relative changes in mass density, (p-p m )/p m for Exp. Bl, 
with respect to the initial density (pi n ), are shown at the same 
times and in the same region as shown in Fig. 6. 



16.0 s 



t = 32.0 s 




t = 64.0 s 



t = 80.0 s 





Fig. 6. Relative changes in the temperature, (T - T m )/T m for 
Exp. Bl, with respect to the initial temperature {T m ), are shown 
at four times in the region 1 < x < 5, < y < 1.6. 



perature for the case of Exp. Bl compared to Exp. Al. This is 
because Joule dissipation becomes stronger as the magnetic field 
strength increases. In the case of Exp. Bl, E\ ~ 800, in the dif- 
fusion region, whereas E2 ~ 500. Even though E2 increases as a 
result of the increase of A(T), the net effect is that more energy 
is available to heat the plasma in the diffusion region. The rel- 
ative changes in the temperature are manifested as a maximum 
in the diffusion region and decrease monotonously in either di- 
rection along the current concentration. The maximum velocity 
of the outflow jets is larger in Exp. Bl than in Exp. Al. This is 
because as the ratio of the magnetic pressure to gas pressure in- 
creases, it leads to higher acceleration of the outflow plasma as 
the magnetic field strength increases. 

Fig. 8 shows the results for the case of Exp. B2, where 
B = 0.033 and p b0 =0.1 (i.e. B = 20GandJV e = 2x 10 10 cnr 3 ). 
In this case, the plasma in the diffusion region is heated to 
?max = 6.4 x 10 5 K. The maximum velocity of the outflow jets, 
Vjet, reaches up to 150 km s _1 , and is 130 km s _1 in the line-of- 
sight. The temperature at the location of vj et is 6 x 10 5 K. Note 
that, for the case with low mass density, the heat conduction term 
will become more and more important as an energy loss in the 
diffusion region compared to the radiative loss term, when the 
temperature of the plasma increases to a higher value. All the 
results for both Exp. A and B are summarized in Table 1. 



Table 1. Summary of the results for both Exp. A and B. v y de- 



notes the line-of-sight component of Vj et . 



1.5 x 10 5 K. For the upflow jet, Vj et reaches 60 km s , directed 
at 2 1 with respect to the line of sight, and Jj et is about 8 x 1 4 K. 
Comparing the results with those in Exp. Al where we use the 
same mass density at the bottom but with a weaker magnetic 
field strength, we find that the plasma is heated to higher tem- 



Quantity 


Exp. Al 


Exp. A2 


Exp. Bl 


Exp. B2 


Tmax (K) 


~1 x 10 4 


2.8 x 10 5 


1.5 x 10 5 


6.4 x 10 3 


Tjet (K) 


~1 x 10 4 


2x 10 5 


8x 10 4 


6x 10 s 


vjet (km s -1 ) 


40 


85 


60 


150 


v y (km s~') 


32 


72 


56 


130 
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t = 32.0 s 




t - 64.0 s 



t = 80.0 s 




SiSsss 



Fig. 8. Relative changes in the temperature, (T - T m )IT m for 
Exp. B2, with respect to the initial temperature {T m ), are shown 
at four times in the region 1 <x<5,0<y< 1.6. 



4. Discussions 

In the present paper, we investigate magnetic reconnection in 
the solar chromosphere, looking for an answer to whether jet- 
like events can originate at lower temperatures. We consider 
magnetic reconnection triggered by newly emerging magnetic 
flux. Several physical environments representing the chromo- 
sphere with different magnetic field strengths and densities are 
discussed. Our results quantitatively show that the temperatures 
and velocities of the bi-directional jets caused by magnetic re- 
connection are strongly dependent on the physical conditions of 
the atmosphere, such as magnetic field strength and mass den- 
sity. In the case of where the magnetic field strength is high 
(Exp. B2) and the density is low, we obtain the highest temper- 
atures and the largest velocities among all our simulations. The 
plasma in the diffusion region is heated up to 6.4 x 10 5 K, and 
the maximum velocity of the outflow jets reaches 150 km s . 
However, in the case of a low magnetic field strength (Exp. Al) 
and high density, the temperature and the velocities of the jets are 
the lowest of all the cases. The temperature is about 10 4 K and 
the maximum velocity of the jets is only 40 km s _1 . Either an in- 
crease in the magnetic field strength or a decrease of the density 
will raise the ratio of the Joule dissipation to the radiative losses. 
That will heighten the heating effect compared to cooling and 
will lead to a high increase in the temperature of the plasma in 
the diffusion region. That will also raise the ratio of the magnetic 
pressure to the gas pressure, resulting in enhanced acceleration 
of the plasma ejected outward from the diffusion region (which 
will lead to the appearance of high velocity jets). 

For the two cases, Exp. A and B, the amounts of the emerg- 
ing flux are not the same, which implies that the stronger the 
magnetic field strength the larger the amount of newly emerging 
magnetic flux. We can also obtain outflows with high velocity 
and high temperature even if a smaller amount of new flux, e.g. 
the same as that in case A, is emerged in case B. Moreover, in 



case B where a larger amount of new flux emerges, the tem- 
perature and velocity of the outflow jets reach their maximum 
at about t = 40 s, not at the end of the flux emergence pro- 
cess. Therefore, the difference in the amount of emerging flux 
between case A and B is not very important in determining dif- 
ferent results in the two cases presented. Instead, the magnetic 
field strength of the initial background magnetic field and the 
initial density structure plays a more significant role. If we do all 
the experiments with constant flux emergence, the amount of the 
newly emerging flux will be very small in case B, compared with 
that of the background field. So that the height of the emerging 
arcade is lower in case B than in case A at any given time, with 
the result that reconnection occurs at different heights in the two 
cases. 

The flow velocities derived here can be compared with 
Roussev et al. (2001a, 2001b) who found blue-shifts of the order 
of 100 km sounder the physical environment of electron den- 
sity 3 x 10 10 cirT 3 and magnetic field 8 G. These authors also 
found that the choice of initial state and the consideration of 
non-equilibrium ionization was crucial to their modelling. Even 
though various initial physical environments were explored in 
their work, the plasma j3 on the current sheet was the same in 
all the cases. Both the plasma p and the thermal energy were 
assumed to be uniform in the y-direction with the reconnection 
happening in the transition region. The maximum jet velocities 
they obtained under different physical environments did differ at 
the begining of the experiments, however, they approached to al- 
most the same value as time went on (see top-left panel in Fig. 9 
in Roussev et al. (2001a)). In our model, the maximum jet veloc- 
ities are different for all the cases. The highest velocity jets are 
obtained under the environment with almost the same density 
as in Roussev et al. (2001a), but stronger magnetic field (20 G). 
When the reconnection occurs in the chromosphere, high veloc- 
ity jets-like features are possible if a larger field strength is used. 
Besides, we do not consider transient ionization. 

In the flux emergence studies of Yokoyama & Shibata 
(1995, 1996) and Moreno-Insertis et al. (2008), reconnection 
occurs at coronal hights where the plasma temperature is high 
and the electron density is low (10 10 cm -3 ), even though new 
flux emerges from beneath the photosphere. Therefore, it is eas- 
ier to obtain jets with high velocities and high temperatures. In 
the Isobe et al. (2008) model, magnetic reconnection between 
newly formed horizontial fields and background vertical fields 
took place at chromospheric heights. The chromospheric elec- 
tron density is much higher (> 10 12 cm -3 ) than that in our model, 
so that the reconnection jets velocity they obtained is very slow, 
about 30 km s _1 . In the Murray et al. (2009) model, oscillatory 
reconnection is driven by the global imbalance of the forces be- 
tween the neighbouring flux systems. At each subsequent recon- 
nection reversal, the system is closer to equilibrium. Once an 
equilibrium state is reached, all reconnection ceases, therefore, 
the maximum temperature and maximum velocity (50 km s~') 
achieved in the outflow jets decreases in each subsequent recon- 
nection phase. All the flux emergence studies mentioned above 
did not include heat conduction and radiative effects. The latter 
will reduce both the temperature and the velocity of the outflow 
jets. 

In comparison to Litvinenko & Chae (2009) model where 
radiative losses were assumed to be of less importance, they 
become significantly important in our model when the den- 
sity and temperature are large. Their study was based on the 
assumption that the inflowing magnetic energy is completely 
converted into thermal and kinetic energies of the outflowing 
plasma. Therefore, the temperature and speed of the outflow jets 
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are larger in their model than those in ours at the same height for 
the same parameter. 

By using our model we would like to give a possible expla- 
nation as to why no brightening is obtained at zero velocity of 
the emission spectral lines in the study by Innes & Toth (1999). 
In their study, the magnetic field strength and mass density in 
the chromosphere was about 6 G and 2.8 x 10~ 10 kg irT 3 , respec- 
tively, which corresponds to (3^ = 1.3 and p^o = 0.8 in our model. 
Our calculations show that in this case the maximum tempera- 
ture in the diffusion region is less than 5 x 10 4 K and therefore, 
no brightening could appear at the zero velocity position of the 
spectral line. 

The radiative loss model that we use at the base is the optical- 
thin approximation. McClymont & Canfield (1983) considered 
optical-depth effects and proposed an optical-thick model for the 
radiative loss function. By comparing the two models, we find 
that the optical-depth effects become important primarily for the 
plasma at temperatures below 3 x 10 4 K, where the contribution 
of hydgron is dominant. This means, in our model, that the radia- 
tive cooling effect is over-evaluated for the cool plasma at tem- 
peratures below 3 x 10 4 K. The temperature of the cool plasma 
should be raised if we take into account the optical-depth effects 
in our model. However, even if we consider optical-depth effects, 
the corrected temperature of the cool plasma will not exceed 
3 x 10 4 K. Otherwise, radiative loss will become too strong and 
the plasma will start to cool again. Considering that the jet-like 
events discussed here correspond to hot jets at transition-region 
temperatures, optical depth effects will have less influnce. 

As mentioned in Section 2, background heating is not con- 
sidered in our model, therefore the initial state is not in thermal- 
equilibrium. We take the case of Exp Al, shown in Fig. 2, as an 
example to examine the change of the background. If no emer- 
gence takes place, the plasma temperature at the top boundary 
will decrease by 27% by the end of 80 s. Because of back- 
ground cooling, downward flow will occur, which reaches about 
40 km s _1 at maximum during 0-80 s. The high velocity region 
is mainly located at high altitude (transition region or above), 
whereas the downward flow around the diffusion region is slower 
(less than 10 km s _1 ). 

Priest & Forbes (1986) studied mechanisms of magnetic re- 
connection analytically. They introduced a new parameter {b) 
and produced a series of analytical solutions. They found that 
the type of reconnection regime and the rate of reconnection de- 
pends sensitively on the parameter b which characterises the in- 
flow conditions. The Petschek (b = 0) solution is just one par- 
ticular member, corresponding to a weak fast-model plasma ex- 
pansion. However, in our model, there is a fast plasma expan- 
sion across the inflow region, and the diffusion region lengthens 
as reconnection occurs. Therefore, the reconnection event here 
resembles more the fast-mode expansion solutions with b > in 
Priest & Forbes (1986). 

In this study, reconnection occurs on the right-hand side of 
the emerging flux arcade. If reconnection occurs on the left- 
hand side of the emerging flux arcade, e.g. the polarity of the 
emerging flux changes, the magnetic field strength on the left- 
hand side will decrease faster as reconnection goes on. This is 
because symmetric boundary conditions are used on the left, 
whereas open boundary is used on the right. As magnetic field 
strength decreases, magnetic reconnection will become weaker. 
Despite of this, there will not be significant differences in the 
results if reconnection occurs on the left-hand side of the emerg- 
ing flux arcade (because the jets reach their maximum velocity 
faster than the decrease of the magnetic field strength). The mag- 
netic field strength of the initial field on the left-hand side of the 



flux emergence region is still very strong, when the jets velocity 
reach their maximum (for instance, at around 40 s in case B). 
If investigating long-term evolution, a larger computational do- 
main, compared with the flux emergence region, is necessary to 
weaken the influence of boundary conditions. 

In the present paper, we mainly concentrate on the hot com- 
ponent of the outflow jets caused by magnetic reconnection, and 
obtain jets with a high temperature and a high velocity under cer- 
tain circumstances. We also produce cool jets in the temperature 
range of (2 + 0.2) x 10 4 K. Taking the case of Exp. B2 as an ex- 
ample. The line-of-sight velocity of the cool upflow jets reaches 
25 km s _1 . This shows one possibility that some cool jets, e.g. 
spicules, and hot jets may originate from the same phenomenon, 
chromospheric magnetic reconnection, and they may just depict 
different parts and/or different stages of the same event. In the 
case of weak magnetic field strength but high mass density, the 
cool jets disappear, which shows that the origination of spicules 
may also have some correlations with the physical environment. 
We relegate further study of the relationship between spicules, 
jet-like events, and magnetic reconnection to future work. 
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